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ABSTRACT 

Recent detections of cosmic microwave background (CMB) anisotropy at half-degree 
angular scales show considerable scatter in the reported amplitude even at similar 
angular resolution. We use Monte Carlo techniques to simulate the current set 
of medium-scale CMB observations, including all relevant aspects of sky coverage 
and measurement technique. The scatter in the reported amplitudes is well within 
the range expected for the standard cold dark matter (CDM) cosmological model, 
and results primarily from the restricted sky coverage of each experiment. Within 
the context of standard CDM current observations of CMB anisotropy support the 
detection of a "Doppler peak" in the CMB power spectrum consistent with baryon 
density 0.01<fib<0.13 (95% confidence) for Hubble constant H = 50 km s -1 Mpc -1 . 
The uncertainties are approximately evenly divided between instrument noise and 
cosmic variance arising from the limited sky coverage. 
Subject headings: cosmic microwave background - nucleosynthesis 
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1 Introduction 



The anisotropy of the cosmic microwave background probes the distribution of matter 
and energy prior to the epoch of structure formation. On large angular scales (0>2°), 
CMB anisotropies reflect perturbations larger than the particle horizon and thus 
probe the primordial density distribution. On angular scales smaller than 2°, causal 
processes become important and modify the primordial distribution in model-specific 
ways. Observations of the CMB at angular scales near 0?5 sample these causal 
processes and offer a powerful test of competing models of structure formation. 

The Cosmic Background Explorer provides a clean detection of CMB 
anisotropy on large angular scales, measuring the amplitude of the fluctuations to 
a significance of 14 standard deviations (Bennett et al. 1996). On smaller angular 
scales, the situation is less clear. A number of groups have reported detections of 
CMB anisotropy on angular scales 0?5 to 1?5 (Figure 1). A large scatter is evident, 
even for experiments at similar angular scales. Does this scatter reflect cosmology 
(e.g. non-Gaussian features in the distribution of CMB anisotropy), instrumentation, 
or simply sample variance from the small sky coverage (typically less than 0.1% of 
the celestial sphere) of an individual measurement? 

Scott, Silk, & White (1995) approach this problem analytically, fitting a 
phenomenological approximation of the CMB power spectrum to recent detections to 
obtain best-fit shape parameters. They conclude that the data do not show excessive 
scatter and obtain best-fit parameters consistent with the existence of a "Doppler 
peak" in the power spectrum at angular scales near 0?5. In this Letter we use Monte 
Carlo techniques to explore further the issue of scatter between different experiments 
and the degree to which instrument noise, calibration error, and sky coverage may be 
expected to affect the results. For specificity, we work in the context of the standard 
cold dark matter cosmology (density — 1, cosmo logical constant A = 0, and 
Hubble constant H = h km s _1 Mpc -1 with h = 0.5), although our results may be 
extended to other choices of model parameters. 



2 Simulations 

We represent the CMB temperature field at angular position (9, 4>) as a sum over 
spherical harmonic amplitudes, T C mb = J2e m a imXim{& \ (f>) ■ The amplitudes a^ m 
are specified by the cosmological model: for a scale-free power spectrum P(k) oc k n 
they are random Gaussian variables with zero mean and variance 

(n 2,_ r 4tt r[Z + (n-l)/2] r[(9-n)/2] 

{a em) — Ue — 5 LJ r[/ + ( 5_ n)/2 ]r[(3 + n)/2]' UJ 

where Q is the ensemble-averaged quadrupole amplitude (Bond & Efstathiou 1987). 
The matter content of the universe will modify the coefficients Ce to increase power 
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at t ~ 200 — 1000 and produce the so-called "Doppler peaks" in the power spectrum. 
In this Letter we use CDM power spectra (Stompor 1994) normalized to the 4- 
year COBE value Q — 18 fiK (Bennett et al. 1996), evaluated over the range 
2f It 500. An experiment with Gaussian beam dispersion <Jb will observe a smoothed 
temperature distribution Tob^^, 0) =J2e m a tmXim{^ ',</>) exp(— ^£(£+l)a$). To reduce 
instrumental artifacts, experimenters rapidly sweep the antenna beam between two 
or more locations on the sky, and combine these to form temperature differences AT 
whose distribution can then be compared to theory. 

A single realization of the CMB sky is specified by the set of amplitudes 
a im- We generate 2000 CMB realizations and use each realization to simulate the 
current set of medium-scale CMB experiments. With a fixed set of amplitudes ae m we 
evaluate T ^ s along the scan patterns and with the beam dispersions of the following 
experiments: ACME/SP94 (Gundersen et al. 1995), ARGO (De Bernardis et al. 
1994), IAB (Piccirillo and Calisse 1993), MAX4 7 UMi (Devlin et al. 1994), MAX4 
cxHer and iDra (Clapp et al. 1994), MSAM 2-beam and 3-beam (Cheng et al. 1994), 
Python II (Ruhl et al. 1995) and Saskatoon '94 (Netterfield et al. 1995). For each 
experiment, we combine the smoothed temperatures T Q b s to form the appropriate set 
of noiseless temperature differences AT (single difference, double difference, and so 
forth). We then add random uncorrelated noise to each set of simulated AT and 
multiply the noisy result by a factor (1 + 5G) where 5G is a Gaussian of mean zero 
and width determined by the calibration uncertainty of each experiment (typically, 
\SG\ ~ 0.1). We remove the mean value from each noisy simulated experiment and 
further remove a linear drift for those groups performing such subtraction. To speed 
processing, we neglect possible foreground emission and correlations from atmospheric 
or detector noise; in cases where the experimenters observe at multiple frequencies, 
we use a single frequency with noise appropriate to the final CMB weighted mean 
used by each group. 

From each sky simulation we thus generate a set of simulated CMB 
experiments, measured at the locations and in a similar manner as the actual 
experiments. Cosmic variance, instrument noise, and calibration uncertainties are 
automatically accounted for. We may then compare the scatter of the simulated 
experiments in each realization to the scatter in the actual data. 

We quantify the amplitude of the CMB anisotropy for each simulated 
experiment using the quadrupole normalization Qflat, the equivalent power for a flat 
power spectrum (Eq. [1] with n = 1) fitted to the simulated temperature differences. 
For each simulated experiment, we derive Qflat by maximizing the likelihood 



where M(Qfl at ) J is the flat-spectrum covariance matrix between pixels a and f3 
computed from the simulations, the index i tracks the different experiments for each 
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realization, and j tracks the 2000 CMB realizations. Note that the subtraction of 
the mean (and sometimes linear drift) from the noisy AT correlates the pixels; we 
account for this by simulating noisy flat-spectrum CMB realizations, computing the 
mean covariance matrix at different values of Qflat, and inverting the (noisy) matrix 
M separately for each Qflat value. 

3 Results 

For a fixed CDM power spectrum, we generate a set of simulated values Q % l &t for the 
ith experiment in the jth realization. We then define a x 2 statistic to quantify the 
scatter in each realization, 

2 / (Qflat ~ (Qflat)) \2 

Xj ~ y " < S Q' "~ ' 

where SQ l is the standard deviation of each simulated Q\ &t evaluated across all 2000 
realizations. The set of 2000 simulated \ 2 values can then be compared to the x 2 °f 
the actual data calculated in a similar fashion (for a recent tabulation see Scott, Silk, 
& White 1995). 

The actual CMB experiments have x 2 = 10-1 f° r 10 degrees of freedom 
when compared to a CDM model with baryon density fixed at tt b = 0.05. 42% 
of the Qb = 0.05 simulations had larger x 2 , placing the current set of CMB 
experiments near the median of the simulated distribution. There is no evidence 
for anomalous scatter in recent detections of CMB anisotropy at half-degree angular 
scales. We repeat the analysis for models with = = 0.2. Aflat "no-baryon" universe 
(Qo = 1, fib — 0, h — 0.5) has \ 2 — 43.0 and is rejected at more than 5 standard 
deviations. CMB anisotropy data clearly support the detection of additional power 
at half-degree angular scales compared to the COBE normalization at large angular 
scales. In the context of CDM with standard recombination, we minimize x 2 to obtain 
a 95% confidence interval for the baryon density 0.01<fi;><0.13. Since we have fixed 
all other cosmological parameters (curvature, cosmological constant, Hubble constant, 
recombination), this range for Qj, can not be regarded as definitive; nevertheless, it is 
reassuring to find the values derived from CMB anisotropy to be in agreement with 
standard Big Bang nucleosynthesis (Walker et al. 1991). 

To what extent do instrument noise, calibration uncertainty, and limited sky 
coverage contribute to the uncertainty in recent medium-scale CMB measurements? 
We repeated the simulations with the same (noiseless) ae m coefficients and varied the 
instrument noise and calibration errors. For noise levels and sky coverage typical 
of recent medium-scale experiments, the combined uncertainties in Qfl a t are in the 
range 6-12 fiK, split almost evenly between cosmic variance and instrument noise; 
calibration uncertainties are generally a minor contribution. Cosmic variance scales 
as the inverse square root of the sky coverage (Scott, Srednicki, & White 1994; 
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Hinshaw, Bennett, & Kogut 1994) and places a lower limit to the accuracy achieved 
at any noise level. Determination of the CMB anisotropy amplitude to 10% accuracy 
at 0?5 angular scales, let alone multi-parameter fits to the power spectrum, will 
require observations with significantly greater sky coverage than is currently achieved. 
Assuming noise levels less than 30 /xK per field of view, a 10% determination of Qflat 
requires observing a 30° x 30° patch on the sky. 
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Figure 1: Recent detections of CMB anisotropy, fitted to the amplitude Qfl at of a 
scale-invariant power spectrum with power equal to the observed power. The grey 
band indicates the 4-year COBE normalization. The solid lines show the predicted 
power spectrum for standard CDM (upper curve) and = (lower line). 
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